Tunneling into d-wave superconductors: Effects of interface spin-orbit coupling 
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Tunneling conductance of a clean normal metal/d-wave superconductor junction is studied by 
using the extended Blonder-Tinkham-Klapwijk formalism. We show that the conductance is signif- 
icantly affected by the interface spin-orbit coupling of the Rashba type, which is inevitably present 
due to the asymmetry of the junction. 

PACS numbers: 74.55.+V, 74.45.+C, 74.20.Rp 



High-temperature cuprate superconductors have re- 
mained at the forefront of experimental and theoretical 
research for more than twenty years. Probing the order 
parameter structure has been the subject of a particu- 
larly large effort. One of the most compelling pieces of 
evidence for the d-wave symmetry of the order parame- 
ter comes from tunneling spectroscopy experiments. The 
conductance spectrum of a junction between a normal 
metal and a high-T c superconductor exhibits strong de- 
pendence on the crystallographic orientation of the inter- 
face, see Refs. [0 andy for a review. Its most prominent 
feature is the zero-bias conductance peaks (ZBCP) that 
can be attributed to the quasiparticle states with zero en- 
ergy bound to the surfaced— Such states, called the An- 
dreev bound states, exist if the quasiparticles experience 
a sign change of the order parameter upon reflection from 
the interface, as was originally pointed out by Hu^ The 
Andreev bound states and the associated low-energy fea- 
tures in the tunneling conductance have also been studied 
for other unconventional superconductors J£~— 

Most of the theoretical studies of the tunneling conduc- 
tance in high-r c superconductors have used the Blonder- 
Tinkham-Klapwijk (BTK) formalism,— extended to the 
(i-wave case. It has been known, however, that going be- 
yond the BTK model produces important qualitative ef- 
fects. For instance, time-reversal symmetry can be spon- 
taneously broken near the interface due to the forma- 
tion of a subdominant order parameter, leading to the 
splitting of the ZBCP even in zero external magnetic 
field) 11 : 12 while in the presence of the interface rough- 
ness, the ZBCP exist for all interface orientations^ Even 
within the BTK framework, the tunneling conductance 
in the (i-wave case turns out to be sensitive to the de- 
tails of the interface barrier, see, e.g., Ref. [l3|, where the 
effects of ferromagnetic and Kondo-like scattering in the 
barrier were considered. 

In this Letter we study the effects of the spin-orbit 
coupling (SOC) localized near an interface between a 
normal metal and a ci-wave superconductor. Due to the 
fact that two sides of the junction have different crystal 
and electronic structure, the interface potential barrier is 
asymmetric, resulting in the SOC of the type originally 
proposed by Rashba in Ref. [lj for semiconductor het- 
erostructures. We neglect disorder as well as the interface 
roughness and calculate the zero-temperatute tunneling 
conductance for different crystalline orientations by gen- 



eralizing the BTK formalism to include the Rashba in- 
terface SOC. Similar model was recently applied in Ref. 
Ha to a normal metal/p-wave superconductor junction. 
Throughout the paper we use the units in which H = 1. 
We consider the tunneling junction shown in Fig. [1] 
The interface is located at x = and is characterized 
microscopically by a potential barrier which we describe 
by the following model: 



U{x) = [U + Um. ■ (<r x k)]S{x). 



(1) 



Here n = x is the unit vector along the interface normal, 
Uq and U\ are the strengths of the spin-independent and 
the Rashba SOC contributions, respectively, <x are the 
Pauli matrices, and k = — iV. The band dispersions are 
assumed to be parabolic, with the same effective masses 
m and the Fermi energies kp/2m on both sides. 

The quasiparticle wave function has four components, 
corresponding to the electron-hole and spin degrees of 
freedom, and can be found from the Bogoliubov-de 
Gennes (BdG) equations^ Assuming two-dimensional 
geometry, the SOC is diagonal in spin, and the BdG 
equations can be decoupled into two independent pairs 
of two-component equations as follows: 



7k*(r) = J5tf(r), 

where a = ± for different spin orientations, 



T~La — 



£ + U a {x) aA{k,r) 
aAt(fc,r) -i-U a (x) 



(2) 



(3) 



| = fe 2 /2m - e F , and U a (x) = (U ~ <rUik v )S(x). The 
gap function is given by A(fc, r) = A(k)9(x), where 9{x) 
is the step function. The spin index in the off-diagonal 
elements of Ti. a amounts to an unimportant phase fac- 
tor and can be dropped when calculating the tunneling 
conductance. We note that the values of the gap can be 
different for the transmitted electron-like and hole-like 
quasiparticles, due to the anisotropy of the d-wave order 
parameter. We have 



A ± = A(fe±) = A cos(26> T 2a), 



(4) 



where a is the angle between the crystalline orientation 
and £-axis, k± = fej? (±cos#,sin#), and A + and A_ are 
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FIG. 1: (Color online) Schematic illustration of the quasi- 
particle reflection and transmission processes at the interface. 
Also shown is the d-wave order parameter profile and the an- 
gle a of the crystalline orientation with respect to a;-axis. 



the effective pair potentials of electron and hole compo- 
nents, respectively. In the spirit of the BTK approach, 
self-conistency of the order parameter is neglected. 

We assume that electrons are injected from the nor- 
mal metal with the excitation energy E > and spin 
a, at an angle 9 from the interface normal. The inci- 
dent electrons are either normally reflected as electrons 
or Andreev reflected as holes.— The momentum parallel 
to the interface is conserved in the tunneling process and 
the solution of Eq. © has the form *(r) = e lfc « y \T/(a;), 
where 



^ N {x) 



1 



ikp cos Ox, 
fc- "T" (*n 



—ikp cos 8 x 



o 



-ikp cos 6 x 



(5) 



on the normal side, and 



Vs(x) 



u+e ^ + 
v+ 

U- 



Akp cos 6 x 



-ikp cos x 



(6) 



on the superconducting side. Here a a and b a are the 
amplitudes of the Andreev and normal reflection, respec- 
tively, and c a and d a are the transmission amplitudes. 
The quasiparticle amplitudes in the superconducting re- 
gion are given by 
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(7) 



where 17 ± = -J ' E 2 — |A±| 2 , with the phase factors e^ ± = 
A±/|A ± |. Note that, according to Eq. (gj), $± = or n. 
All the reflection and transmission amplitudes in Eqs. 
([5]) and © can be found from the boundary conditions 
that follow from Eq. ([1]): 

*s(0 + ) = *at(0-), 

*' s (0+) - ^(0-) = 2m(U - aUik F sin9)^ N (0-). 
In particular, for the reflection amplitudes we obtain: 

4r+e _i *+ 



a a {E,9) = 

ba{E,9) = 

where uj± = 1 ± r + r_ 
depending on 9 and a), 



r+ = 



z a = 



V± = |A±| 
u± E + tt± 
Zq — aZ\ sin 9 
cos 9 



(2 + Z 2 )u;_ + 2uj+ ' 

{-2iZ a - Z2)u- 

(2 + Z 2 )lu- + 2uj+' 



e-q^ 



(8) 



(e^ = ±1 



I Ah 



Zn = 



2mU t 



o 



&F 



Zi = 2mC/i. 



The dimensionless parameters Zq and Zi characterize the 
strengths of the purely potential and SO scattering, re- 
spectively. 

Using the BTK formalism— the normalized dif- 
ferential tunneling conductance is given by G{E) — 
G S (E)/G N , where 



G S (E) 



E 



7T/2 



71-/2 



dO cos 9 G a (E, 



(9) 



with the angle and spin resolved conductance given by 

I 2 
-| 2 )/4 



G a {E,9) = l + \a a (E,9)\ 2 -\b a (E,9) 

i + |r + | 2 + ^(i-|r + r_ 



|i + z2(i-r+r_e^)/4| 2 



and 



Gn — 



tt/2 



d9 cos 9 



1 



tt/2 



1 + Zl/A 



(10) 



(11) 



is the conductance for a normal metal/normal metal 
junction with the interface potential given by Eq. (JlJ. 
We can see that the tunneling conductance depends on 
the incident spin orientation: G+(E,8) ^ G-(E,9). In 
the absence of the interface SOC, i.e. at Z\ = 0, Eq. (fT0| 
reduces to the known results [see, e.g., Eq. (23) of Ref. 

H- 

The effect of the Andreev bound states is most pro- 
nounced in the low-transparency limit, i.e. when the in- 
terface barrier is so high that one can put Z a — > oo. In 
this case, G^ becomes small, but G a (E,9) remains inde- 
pendent of the barrier height if r + r_e J< ^ = 1. The last 



equation can be written in the form A + /(E — fl+) — 
A_/(E + fi_), which has a zero-energy solution at 
a = tt/4 for all incident angles due to the fact that 
A_ = — A + (Ref. [3|). Therefore, the normalized tun- 
neling conductance at E = diverges, giving rise to a 
sharp ZBCP. In contrast, at a — there are no zero- 
energy bound states near the interface and the tunneling 
probes the density of the quasiparticle states N(E) in the 
bulk. In the d-wave case, at low energies the main contri- 
bution to the density of states comes from the vicinity of 
the gap nodes, yielding N(E) ~ E (Ref. [H). This leads 
to a strong suppression of the tunneling conductance at 
low bias. We shall see below that the last conclusion 
surprizingly changes if the interface SOC is taken into 
account. 

We use expressions © and pT|) to calculate the nor- 
malized tunneling conductance for different values of the 
interface SOC at a — and a = ir/4. For the po- 
tential barrier height we consider three cases: Zq = 
(high transparency), Zq — 1 (medium transparency), and 
Zq = 5 (low transparency). The conductance is plot- 
ted as a function of the dimensionless excitation energy 
E/A . 

Figures 2 and 3 show the tunneling conductance for 
a = and a = 7r/4, respectively, in the high and medium 
transparency cases. We note that if Zq = Z\ = 0, then 
there is only Andreev reflection and G{E) is nearly in- 
dependent on a, monotonically decreasing from 2 as £ 
increases. This changes in the presence of the interface 
SOC, when a maximum appears in G{E) in the subgap 
region at a — 0, as shown in the top panel of Fig. 2. 

Another significant feature of our results is that in the 
high-transparency case the effects of the SOC on the zero- 
bias conductance are opposite for the two interface ori- 
entations: At a — G(0) is suppressed as Z\ increases, 
while at a = ir/4 G(0) is enhanced, as shown in the top 
panels of Figs. 2 and 3, respectively. To explain this we 
use the analytical expression for the zero-bias conduc- 
tance at Zq = 0: 



G(0) = 



(x\ + l)(xf — 3) arctan^i -I- {x\ + 3)xx x\ 
(x\ + 1) arctanx2 — X2 x\ 



where x\ = ^/ Zf{l + e^/A- 1 and x 2 = ^Zf/4-1. 
At Z x = 0, we find G(0) = 2. At a = (a = tt/4), we 
have e 1 ^ = 1 (e'^ = —1), and the above expression is a 
decreasing (increasing) function of Z\ . 

In the medium-transparency case the interface SOC 
produces some enhancement of the tunneling conductiv- 
ity, which is more pronounced in the a = case. 

Let us now discuss the case of a low transparency bar- 
rier, which is shown in Fig. 4. At a — tt/4, the interface 
SOC suppresses the height of the ZBCP. At a = 0, there 
is no ZBCP, but the tunneling conductance at a small 
bias is enhanced by the SOC, reaching values of the or- 
der of 1 (in contrast, the peak at E = A is suppressed). 
To explain this enhancement we note that for Zx smaller 
than Z we have Z a ~ Z Q > 1. Then G s (0) - Z ~ 4 , 
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FIG. 2: The dimensionless tunneling conductance for a = 
and different strengths of the interface SOC, at Zq = (top 
panel) and Zo = 1 (bottom panel). 




FIG. 3: The dimensionless tunneling conductance for a — n/4 
and different strengths of the interface SOC, at Zo — (top 
panel) and Zo = 1 (bottom panel). 
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while Gjv ~ Z Q 2 , therefore, (7(0) - Z 2 -> 0, sim- 
ilarly to the zero SOC case, see, e.g., Ref. |4|. How- 
ever, if Z\ > Zq, then for some incident angles one has 
\Z a \ <C 1. It is the contribution from those angles that 
dominates the integrals in Eqs. (O and (jlip . giving rise 
to G s {0),G N ~ Zf 1 , therefore (7(0) ~ 1. 



To summarize, we have calculated the tunneling con- 
ductance of a junction between a normal metal and a 
d-wave superconductor. Unlike the previous works, we 
take into account the SOC localized near the interface, 
which requires a modification of the BTK formalism. We 
have shown that the interface SOC gives rise to sev- 
eral qualitative changes in the tunneling spectra. The 
most prominent changes are as follows: In the case of 
a high-transparency junction, the normalized tunneling 
conductance at a small bias E <C Ao is suppressed by the 
SOC for a — and enhanced for a — 7r/4. In the low- 
transparency junction, the trends are reversed, in partic- 
ular, the zero-bias conductance at a = is enhanced by 
the SOC. 



FIG. 4: The dimensionless tunneling conductance in the low 
transparency case, for different strengths of the interface SOC, 
at a = (top panel) and a = tt/4 (bottom panel). 
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